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Disconnected diagrams give crucial contributions to the physics of flavor singlet hadrons and to 
scalar form factors of non-singlet hadrons. Naive lattice calculation of the disconnected diagrams, 
however, requires a huge number of fermion matrix inversions and hence a prohibitively large 
computational cost. In this article, we present recent studies of the flavor-singelt meson spec- 
trum and nucleon strange quark content using the all-to-all propagator to calculate the relevant 
disconnected diagrams. 
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Lattice studies ofhadron physics with disconnected quark loops 



1. Simulation method 

Disconnected diagrams give crucial contributions to interesting physics including the flavor- 
singlet meson spectrum and nucleon strange quark content. Interesting nature of the former, such 
as the famous U(l) problem and the ideal mixing of vector mesons, has to be confirmed from first 
principles, and the latter is an important parameter in direct experimental searches for the dark 
matter. In this article, we present studies by the JLQCD Collaboration [jl|, ^, |3]] that use a recently 
proposed method, namely the all-to-all quark propagator 

Our gauge configurations of three-flavor QCD are generated at a lattice spacing a ~ 0. 1 1 fm 
using the overlap quark action, which exactly preserves chiral symmetry on the lattice. We simulate 
four values of degenerate up and down quark masses mi, which cover a region of the pion mass 
300<M^[GeV] <550 MeV, and two strange quark masses m s ^m s , p hys and (5/4)m SjP hy S withm ijP hy S 
the physical mass. The flavor-singlet meson spectrum is calculated on a spatial lattice volume 
(1.8 fm) 3 , and we also simulate a larger volume (2.7 fm) 3 for the nucleon strange quark content to 
avoid finite volume effects. We refer readers to Ref. [Q] for further details of our gauge ensembles. 

It is conventional to calculate connected hadron correlators by using the so-called point-to-all 
quark propagator Sp{x' ,x), which flows from a fixed lattice site x to any site x' and is obtained by 
solving 

Y,D(y,x')S F (x',x) = 8 y , x , (1.1) 
y 

where D is the Dirac operator. One can naively evaluate disconnected quark loops Sp(x,x) by 
solving Ql.l[ ) for each lattice site x, but it requires prohibitively large computational cost. Instead, 
we may construct the all-to-all quark propagator, which contains propagations from any lattice site 
to any site, in an effective way [Q]. We determine low-lying eigenvalues {X{\ and corresponding 
eigenvectors of D, and the contribution of these low-modes to the all-to-all quark propagator 
is calculated exactly as SF t i ow (x,y) = Y!kL\ ~t u k{ x )u\{y)- The number of the low-modes N e is set 
to 160 (240) on the smaller (larger) lattice. We confirm that low-energy observables, such as 
disconnected hadron correlators, are well dominated by the low-mode contribution with this choice 
of N e . The small contribution from the higher-modes can be estimated stochastically by the so- 
called noise method, which is not computationally intensive. This all-to-all propagator is used to 
calculate disconnected correlators of the flavor-singlet mesons and nucleon. We refer to Refs. [|l|, 
^[ |6|] for technical details. 

2. Flavor-singlet meson spectrum 

In the left panel of Fig. [j], we plot vector meson correlators Gy^bi^t) = (I/A^L^Vq^ + 
At)V^(t)) of the light (V/ = \fllykl) and strange vector mesons (Vs = s%s) as well as those of 
the flavor-octet (V$ = y/ 2/3(1^1 — sYks)) and singlet mesons (Vq = (2%/ + syts) /\/3). The off- 
diagonal correlators Gy {/ v /j of V/ and V s are about two orders of magnitude smaller than the 
diagonal ones G v sii 7SS \, while there is no such large hierarchy between Vq and Vg. We can express 
CO and <p meson states as eigenvectors of the 2x2 correlator matrices Gy.ab (a,b€: {l,s} or {8,0}) 

= 1.00(1) V s - 0.04(9) V, = 0.84(5) V 8 - 0.55(7) V (21) 
co = 0.04(9)V S + 1.00(1)V/ = 0.55(7) V 8 + 0.84(5) V ' 
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Figure 1: Left panel: vector meson correlators Gy „/,. Right panel: chiral extrapolation of M^, . The horizon- 
tal line represents the bare light quark mass in lattice units, and m S ph ys is 0.08 1 in this normalization. 

These observations imply an almost ideal mixing of the vector mesons, which is expected from the 
experimental spectrum. The right panel of Fig. |l| shows a linear chiral extrapolation of the mass of 
the <p meson state in (2.1). At the physical quark masses, we obtain = 1036(12) stat (— 91) sys MeV 
and M m = 854(24) stat (— 92) sys MeV, where the systematic error is estimated by adding higher order 
terms into the extrapolation form. 

A similar analysis for the PS mesons suggests that 77 and v\' mesons are significant mixture of 
the light (Pi) and strange mesons (P s ) in contrast to the vector mesons: 



17 =0.96(1)^-0.28(3)^, T]' = 0.28(3)^ -0.96(1) P s . 



(2.2) 



A linear chiral extrapolation yields M v = 620(47) sta t MeV and = 789(128) sta t MeV. The 
obtained spectrum of the vector and pseudo-scalar mesons is in good agreement with experiment. 



3. Nucleon strange quark content 

We extract the bare value of the strange quark content (A r |^|/Y)b a re from the ratio of the three- 
and two-point functions 



R(At,At') 



Y t (0if(t + At + At')\{ss(t + At)} hm \0 N (t)) 



Af,Af'- 



-> (N\ss\N) 



bare' 



(3.1) 



A clear non-zero signal shown in the left panel of Fig. |2] is obtained by improving our simulation 
set up [ph: namely, we use the all-to-all propagator to precisely calculate the disconnected quark 
loop, and also use smeared nucleon source (@n) and sink operators (^v) in order to reduce excited 
state contamination of the three- and two-point functions. 

In our previous study in two-flavor QCD [0], we observe that the heavy baryon chiral ex- 
pansion of {N\ss\N)bare shows poor convergence at similar quark masses. We therefore employ 
a simple linear extrapolation to the physical point and also test a constant fit to estimate the sys- 
tematic uncertainty. Those fits are shown in the right panel of Fig. ||. The extrapolated value is 
converted to a renormalization group invariant parameter 



m s (N\ss\N) 

Mm 



= 0.013(12)^(16), 



sys • 



(3.2) 
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Figure 2: Left panel: ratio R(At , Af') with Af+Af'= 13a as a function of At. Right panel: chiral extrapolation 
of (A^|ii|A^)bare- Solid and dashed lines show linear and constant fits at m s —m s iP h ys . 



We note that chiral symmetry, which is preserved in our simulation, is crucial to avoid a possibly 
large contamination due to the operator mixing with the light quark content (N\iiu + dd\N\ 
We also emphasize that (3.2) is in good agreement with our previous estimate in two-flavor QCD 
[^] and our indirect determination through the Feynman-Hellmann theorem [{J: all our studies 
consistently favor the small strange quark content fj s ssO.01 -0.03. 



4. Conclusion 



In this article, we demonstrated the feasibility of a quantitative study of the flavor-singlet 
spectrum and nucleon strange quark content by using the all-to-all propagator. It is interesting to 
extend this study to other observables, such as the nucleon spin fraction carried by sea quarks. 

We also note that the pion scalar form factor has a significant contribution from the discon- 
nected diagram [Q]. Although disconnected contributions are generally (much) smaller than con- 
nected ones, they can not be naively ignored in future precise determinations of hadron observables. 

Numerical simulations are performed on Hitachi SRI 1000 and IBM System Blue Gene Solu- 
tion at High Energy Accelerator Research Organization (KEK) under a support of its Large Scale 
Simulation Program (No. 09/10-09). This work is supported in part by the Grant-in-Aid of the 
Japanese Ministry of Education, Culture, Sports, Science and Technology (No. 20105005 and 
21684013). 
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